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' Abstract 

Recent analytical and numerical results for the three-loop polarization function allow to 
present a phenomenological analysis of the cross section for massive quark production 
in electron positron annihilation to order a^. Numerical predictions based on fixed order 
perturbation theory are presented for charm and bottom production above 5 and 11.5 GeV, 
respectively. The contribution from these energy regions to a(M^), the running QED 
coupling constant at scale Mz, are given. The dominant terms close to threshold, i.e. in 
an expansion for small quark velocity /3, are presented. 
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1 Introduction 



The total cross section for e^e~ annihilation into hadrons, ahad, constitutes one of the most 
basic quantities of hadronic physics. It can be determined experimentally and calculated 
theoretically with high precision. It allows for a fundamental test of QCD and for a precise 
determination of its parameters, the strong coupling constant, and the quark masses. In 
addition, it provides the decisive input for an evaluation of the running QED coupling at 
high energies and for the hadronic contribution to the lepton anomalous magnetic moment. 
Perturbative QCD is expected to provide reliable predictions in the continuum, i.e. one 
or two GeV above the respective quark threshold and the respective resonance region. 
Calculations in the massless limit have been performed several years ago in 0(a;^) and 
0{al) [0 which we also refer to as NLO and NNLO (for a review see [01). The effect of 
nonvanishing quark masses, Mq, has been taken into account during the last years with 
the help of various quite different approaches: a large number of terms has been calculated 
in an expansion in Mq/s [|[ ||, where -/s is the center of mass energy, and a subset of 
the diagrams has been evaluated analytically p, . Alternatively, the real and imaginary 
part of the polarization function n(g^) has been obtained by deriving analytical results 
for the expansions around = 0, for the limit Mq/q^ <C 1 and around = 4:Mq and 
by using the analyticity of 11 (g^) to reconstruct the full function numerically ||^. The 
polarization function to order for massive quarks is therefore under full control. 

The previous papers were devoted to the technical aspects of the calculation and to 
systematic tests and cross checks. The present paper will be devoted to a compilation of 
the results in a simple and coherent form and to various phenomenological applications. 
It contains, in addition, the contribution from the double bubble diagram with a massive 
quark in the internal and the external fermion loop. In the phenomenological applications 
the dominant terms of order in the massless approximation plus Mq/s terms ^ will 
be included. This approach allows for a smooth interpolation between the high energy 
region where the formulae are accurate to NNLO order and the region closer to threshold 
where the results are valid to NLO order only. 

The paper is organized as follows: a comprehensive account of all NLO results is 
presented in Section ^ for the sample case of the charm cross section. In Section ^ 
predictions for charm, bottom and top quark cross sections will be given. The sensitivity 
of the results towards a variation of the input parameters and the renormalization scale 
is investigated. In view of their stability in the regions away from the resonances the 
results for the cross sections can be used to predict the contributions of the charm and 
bottom continuum to the running of the QED coupling. A detailed study of this effect is 
performed in Section ^. The NLO perturbative results expanded for small velocities are 
an essential input for a calculation of the cross section very close to threshold, i.e. for 
energies comparable to or smaller than the Rydberg energy. In this region a resummation 
of leading and subleading terms of order nas/ /3, where [3 = ^1 — 4Mq / s is the velocity 
of the produced quarks, is required. Although we do not perform this resummation 
of singular terms in this paper, the essential ingredients from perturbation theory are 
presented in Section ^. Section |^ contains the summary and conclusion. 
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2 Predictions of order 



In a first step the theoretical results shall be recalled which are required for the complete 
prediction of order valid for energies sufficiently above the heavy quark threshold. The 
crucial ingredient in the present approach is the existence of a hierarchy of the quark 
masses. To be specific, let us consider the region above the charm and below the bottom 
threshold — the generalization to the other cases of interest being obvious. The energy 
is chosen sufficiently large, say above 5 GeV, to avoid the complications in the regime 
very close to the cc threshold (see Section |^). The u, d and s quark masses are neglected. 
Virtual bottom effects are treated through an expansion in s/(4M^). This approximation 
is adequate in the full energy region under consideration, even for s/(4M^) ^ 1 We 
are mainly interested in the region where and Mc are of comparable magnitude. It 
is thus convenient to identify the quark mass with the pole mass, a convention adopted 
throughout this paper. 

In order the contributions from different quark species to the vector current cor- 
relator can be distinguished and thus added incoherently — only "non-singlet terms" 
are present in this order. The singlet contribution, which starts in order a^, has been 
calculated for massless quarks and is small 0. In the following we shall first recall the 
contributions arising from the electromagnetic current coupled to the light u, d and s 
quarks, and subsequently the charm contribution — the main subject of this work. We 
would like to stress, that the formulae, with the obvious replacements, are equally well 
applicable for bb or tt production above their respective thresholds (see Section H). 

1. The sum of the absorptive parts of one-, two- and three- loop diagrams with massless 
degrees of freedom (quarks or gluons) is given by |[l|]: 



where C(3) ~ 1.2020569 and = — 1 is the number of massless quarks. The MS 
coupling af"* is to be evaluated at the scale /i^. Anticipating our strategy to include 
the dominant terms in the high energy region, the evolution of the strong coupling is 
governed by the three- loop beta function with nf = ni + 1 active fiavours, where Uf = 4 
is chosen for the sample case of charm quark production. 

2. Charm quarks can be produced through the splitting of gluons, which in turn are 
radiated off u, d or s quarks. The analytic result for this cross section can be found in |0, 
the corresponding virtual corrections to light quark pair production were obtained earlier 
in [jlO| in the on-shell renormalization scheme. The sum gives rise to the following "double 
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Figure 1: Fermionic double bubble diagrams with generic masses mi and m2. 



bubble" contribution (see Fig. |l| with mi = 0, 1712 = Mc) 



TCRfJ = 3 E Q'^l(p'^iO,Mls) + p\0,Mls) + -ln^] 

i=u,d,s \ ^ / 



(2) 



The functions p^{0,M^,s) and p^(0, M^, s) are given in 0. The combination 
on the right hand side of Eq. is well approximated by the leading terms in 
the high energy limit. This is demonstrated in Fig. ^, where the combination 
(^p^(0, M^, s) + p^(0, M^, s) + ^ In is shown, together with the leading terms of the 
high energy approximation (/i^ = s) 
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Note that the absence of terms could be inferred from general renormalization group 
considerations 0. 

3. Double bubble diagrams with external u, d, s (or c) and internal b quarks decouple in the 
limit s/(4M^) ^ 1. For vanishing external quark mass an analytic result is available |T0| , 



closely related to p (0, M^, s) given above. It is well approximated by the leading term 

(2) ^ ^ S 



in the s/M^ expansion |[TT| , |T2[ — even up to s = 4M^- 



i=u,d,s b 



8 , 176 
In — H 

135 s 675 



(4) 



These terms are numerically small. 

The combination of all terms proportional to J2i=u,d,s Ql thus reads 

Rudsis) = Rii,u+(^^^] TC^(i?g(s)+54?(s)). (5) 
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Figure 2: The function (p^(0, M^, s) + p^(0, M^, s) + ^ In ^) as described 
in the text. Sohd hne: exact result; dash-dotted line: constant and quadratic 
terms only; dashed line: including terms up to order M^/s"^ corrections 
(from [H). The scale ^'^ = s has been adopted. 




Ruds is separately renormalization group invariant and approaches -Riightln^^n^ in the limit 
M2 < s < Ml 

Let us now proceed to the contributions arising from charm quarks coupled to the 
electromagnetic current. They will be cast into the form 

= Ql [r^^ + C, i?« + ' Ri?) 1 . (6) 

The lowest order terms are well known 1113 and read 



= 3/3^, = 3p^'\ (7) 



where 

2 U2{p) + Li2(p^) + \np (2 ln(l - p) + ln(l + p 



^ 2 



- /3 (3 - r ) (2 ln(l - p) + ln(l + p) 

(1-/5) (33-39/3-17/32 + 7/33) ^ 3/3 (5-3/32) 

mp H 
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with 



P 



1-/? 
l + (3 



(3 = ^l-AM^/s 



(9) 



and Lin,(p) is the polylogarithmic function. In the hmit /5 — p^^^ behaves as follows: 



P 



(1) ^ 
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C(2)-6/3 + 3C(2)/32 + 0(/53 



(10) 



with C(2) = 7r^/6. In order a variety of diagrams has to be considered. 



4. The essential ingredients for an evaluation of the double bubble diagram with two 
charm quark loops (see Fig. |l| with mi = 1712 = Mc) can be found in 0. The virtual 
corrections to the cc vertex contribute for ^/s > 2Mc and are known analytically . The 
final state with four charm (anti-) quarks is strongly suppressed close to its threshold at 
4Mc ^ 7 GeV. The rate is given in terms of a two dimensional integral to be evaluated 
numerically. The combined contribution is thus written as 

TC,R^^ = 3 (^pW,M,^.) + ^p^(M,^M,^.) + lln^^p«) (11) 



where 

p"" (Ml Mis) 



3 + 10/32 - 5/3^ 



In^p + 



-3 + 40/32 + 16/3^ - 15/3^ , 2 



+ 
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-18 + 234/32 + 167/3^ - 118^6 -3 - 10/32 + 5/3^ 



I8/32 

9 + 510/32 - 118/3^ 



and 

P«(M2,M2,.) 

with 



3 JiM^/s 
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+ /3(-27 + 5/32)C(2) 
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C(2) ) Inp 
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4M2 /s 2; 



1 + 



2M2 



sz 



I 4M2 , , 

1 '-Hy,z) 
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(13) 
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c ? ■ 



f2 , 



vanishes 



Note that p^{M^,M^,s) vanishes for s 16M^. The function p^{M- 
for s 4M^ and p^^^ approaches the constant 37r^/4 in the same hmit. Both are zero 
below 4M^. For small (3 one obtains: 



(f-4C(2)),a + 



245 8 



54 



+ ^C(2) ) /3' + 



(16) 



R'^ is shown in Fig. ^as a function of M^/s in the range from to 1/4. The contribution 



Mls) 



+ 



from four particle production p {M^, M^, s) and the virtual correction p {M^ 
ip(^) In — are displayed separately as dashed and dotted lines, respectively, their sum is 



shown as a solid line. In the high energy limit the sum approaches (for p"^ 
value: 



s) a constant 



1 A/f2 

p^'iMl Ml s) + p"{Ml Ml s) + ipW In 



C(3) -'i-^oim. (17) 



Fig. 1^ demonstrates that for energies far above the four particle threshold, i.e. for 
^/s 3> 4Mc, real and virtual contributions cancel to a large extent. For smaller ener- 
gies, however, the (negative) virtual corrections become increasingly more important as 
the energy decreases down to the two particle threshold. 

The contributions from charm quarks coupled to the external current with internal 
massless quark or gluon lines are significantly more important. Their treatment is the 
main subject of this paper. Very close to threshold the Coulomb singularity has to be 
incorporated and the definition of the coupling has to be scrutinized. However, in a first 
step, the energy region will be considered where mass terms are important but Coulomb 
resummation is not yet required. 

5. Let us start with double bubble diagrams with light internal quark loops (see Fig. ^ 
with mi = Mc,m2 = 0). In the previous cases, Eqs. (0) and ( pH]) with massive internal 
quark loops, the rates for real and virtual radiation could be given separately and no mass 
singularity was present. This differs from the case with vanishing internal quark mass: 
quadratic and linear mass logarithms arise in the individual cuts which can be cancelled 
by combining real and virtual emission and by adopting the MS definition of the coupling 
constant. The analytical result has been obtained in (see also fl^). For completeness 
we recall the result for light quark species: 



ne 



TCf r9) 



cq 




where p^^^ is given in Eq. (^ and 
^p, _ _ (3-, 9^) 



s 3 



In 4 



• I Li3(p) - 2 Li3(l - p) - 3 Li3(p2) - 4 Li3(Y^) - 5 Li3(l - 
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Figure 3: Second order contributions from double bubble diagrams with 
massive internal and external quarks of the same mass as a function of M^/s. 
Dashed line: the real contribution p^{M^, M^, s); dotted line: the virtual 
correction p^{M^, M^, s) + ^p^^^ In ^ for p"^ = s; sohd line: the sum of both. 



+ Li2(p) ln( 



-)+2Li2(/)ln( 



2/32 



- In (^^) 36 ln2 Inp - 44 In^p + 49 Inp In — —) + In^ ( 



)+2C(2) lnp-ln( 



1-/9' 



2 -I 



4/9 



- Inp ln/3 



36 ln2 + 21 lnp+16 ln/3-22 - p^) 



+ ^ I (15 - 6/3^ - P') {U,{p) + U,{p')) + 3 (7 - 22/32 + 7/3^) Li^b) 
- (1 -/3) (51 -45/3- 27/32 + 5/3=^)C(2) 

+ (-9 + 33/3-9/3^-15/33 + 4/3^) In^p 



+ 



/3 

(33 + 22/3^-7/3^) ln2 - 10 (3 - /32) (1 + /?2) In/? 
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Figure 4: Gluonic double bubble diagrams where the external photon is 
coupled to a heavy quark and the emitted gluon splits into a gluon loop. The 
diagrams for the ghost particles which also have to be taken into account, are 
not depicted. 



-(15-22/5^ + 3/3^) In(^) 



+ 2/3(3-/32) ln( 



4 (1-/3^ 



4/3 

In /3 - 3 In ( 



Inp 

1-/3^ - 
4/3 



+ 



237 - 96 /3 + 62 /32 + 32 /33 - 59 (3^ 
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Inp- 16/3(3 -/32) ln( 



2/3(39- 17 /32) ln(- 



1-/3^^ /3 (75 - 29/32) 



For small velocities 5*^2) jg given by: 



1 + /3 



3C(2) In 



/5 



I + 8 In 2) /3 + 2C(2) (in ^-2) /32 + 0{P') 



(19) 



(20) 



6. Diagrams with massive quarks and purely gluonic internal lines have been evaluated 
in P] through a combination of analytical and numerical methods. The decomposition of 
the result according to the colour structure will be important for the discussion in Section ^ 
below. Terms proportional to with a threshold singularity proportional to {TxasY' / 1^ 
present in abelian and nonabelian theories as well, whereas terms proportional to CfCa 
with a logarithmic threshold singularity are characteristic for the nonabelian structure of 
the theory, with a behaviour similar to the CpTni term. This leads to the decomposition 



CpR^ 



(2) 



f211 



(2) 

where 5-Rc6 is the contribution with an internal b quark loop obtained in analogy to 
Eq. (§. 
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The following approximations have been derived in 

.2) (1-/52)^371^ _ ^2619 ^, 2061 81 /, ^^X , 

/ /\/f2\ 3/2 

_198 (/34-2/32 



100^3/2(1 - p) (2.08 - 1.57p + 0.405p2) , (22) 



/?(2) _ d(2) 



«=4 ^ 32 ^ 32 



-18 (^i^ j - 2/3^ 

+50p^/2(l - p) (l.41 - 1.24p + 0.96p2) , (23) 
where p = (1 — /3)/(l + /3). R^^^\^=4: is the contribution from gluonic double bubble 



diagrams (see Fig. |) for the special choice of the gauge parameter ^ = 4 and reads |jT5 



R(2) = fii In + p(i) _ ^ 5(2) , (24) 

with p(^) and 5*-^'' given in Eqs. (H) and ([T9|), respectively. The first lines of Eqs. (^) and 
(p3|) consist of the exactly known high energy and threshold contributions, the second 
and third lines represent a numerically small reminder. 

3 Cross section for the heavy quark production 

The collection of the results presented in the previous section provides all tools necessary 
for a complete description of the cross section in NLO, including charm, bottom and top 
mass terms. This allows for the prediction of the charm, bottom and top cross sections 
in the regions where quark masses cannot be neglected but where the resummation of 
Coulomb terms, characteristic for the regime very close to threshold, is not yet necessary, 
see the discussion below. As stated above, the terms proportional to Ql and J2i=u,d,s Ql 
are invariant under renormalization group transformations separately, and only terms 
proportional to Ql will be considered in the following. Mutatis mutandis the same for- 
mulae are applicable to bottom quark production. For top quarks only the piece induced 
by the electromagnetic current will be considered, the axial part has been calculated 
in H, m, 0. 



Fixed order perturbation theory is inapplicable very close to the production threshold, 
in the region where as/ (3 is of order one or larger. In this region terms proportional 
to (as/P)"' have to be resummed. However, it will be demonstrated in Section |^ that 
the first three terms in the perturbative expansion provide a good description down to 
fairly small values of f3. Specifically, the relative deviation amounts to 0.9/1.7/3.3% for 
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CpT^ots/ (3 = 2/2.5/7r, respectively. These values lie well within the radius of convergence 
Cpnag/P < 2tt of the resummed series. Taking the requirement Cpnas/P < 2 as a 
guiding principle and incorporating the running of the coupling constant one would admit 
the strictly perturbative, fixed order treatment down to energy values which are 1 GeV 
above the nominal threshold for bottom quarks and even less for charm quarks. However, 
since the perturbative treatment can only be applied beyond the highest cc and hh bound 
states, we take 5 GeV for charm and 11.5 GeV for bottom quarks as lowest center of mass 
energy values. For top, on the other hand, the limit Cpiras/ (3 < 2 corresponds to energies 
about 12 GeV above 2Mf. 

The compensation between phase space suppression and Coulomb enhancement leads 
to a fairly flat energy dependence of R{s) even relatively close to threshold (see Fig. H). 
In [H, 1^ it has been demonstrated for 0{as) and 0{a1), respectively, that this behaviour 
is well approximated by the first terms in the large momentum expansion. Assuming that 





Figure 5: The functions Re, Rb and Rt in NLO plus dominant NNLO terms 
versus y/s for three different scales, /i^ = Mq (dashed), /i^ = (2Mq)^ (solid) 
and fi"^ = s (dotted curves). For comparison, also shown are the Born (wide 
dots) and 0{as) results (/x^ = (2Mq)^, dash-dotted). 
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(a) (b) 

Figure 6: Variation of q;^^-'(M|) (a) and the quark masses (b). In (a) the sohd, 
dashed and dotted curves correspond to ai^\Ml) = 0.118,0.115 and 0.121, 
respectively. In Figs, (b) for the solid curves (Me,Mfe,Mt) = (1.6, 4.7, 175) GeV 
is chosen. The dashed curves correspond to the upper ((1.8,5.0,180) GeV) 
and the dotted curves to the lower limits ((1.4, 4.4, 170) GeV). The scale /x^ = 
{2MqY has been adopted. 



the same line of reasoning is applicable also in order a^, the dominant NNLO corrections 
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Figure 7: R{s) plotted against y/s. The scale jj^ = s has been adopted. The 
dashed curves correspond to the values Mc = 1.8 GeV, Mb = 5.0 GeV and 
as(M|) = 0.121, whereas for the sohd curves = 1.4 GeV, Mb = 4.4 GeV 
and as(M|) = 0.115 is used. The dotted lines show a recent compilation of the 
available experimental data. The central curves correspond to the mean val- 
ues, upper and lower curves indicate the combined statistical and systematical 
errors. 



can be incorporated B 0, IT 



R'^^l, = sj - 6.637+ 17.296 In— + 7.563 In^^ 
[ s s 

+ nf 1^-1.200 - 2.088 In ^ -0.9171n^^^ 

+ n^f -0.005 + 0.038 In ^ + 0.028 In^ ^ + ^ 
V s sis 



347.168- 378.000 In 



M2 
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- 9.000 In^ At + 974.250 In — - 114.000 In At In — + 213.750 In^ — 

M2 s s s 

+ Uf ( -67.619 + 17.333 In Ar + 2.000 In^ A - 82.1671n — 

2 2 2\ 

+ 4.000 In A In — - 17.000 In^ ^ 

s s J 

+ n) ( 1.218 + 1.444 In y + 0.333 In^ ^ 



(25) 



Inclusion of these terms will lead to the correct NNLO predictions for larger energies, say 
above 7 or 8 GeV, allowing at the same time for a smooth transition to NLO accuracy 
for lower energies. The singlet terms proportional to {^QiY are small and have been 
neglected in Eq. (pSf). In total one thus finds: 

fl, . Q?U + ^C.fl«H-(^)'flf)+(^£Myflgi (26) 

where R^^ns denotes the non-singlet contribution at 0{al). In Fig. || the predictions 
are shown for the charm, bottom and top cross sections in the energy regions discussed 
above. For the value of the coupling af\M'^) = 0.118 has been adopted and the on- 
shell quark masses have been chosen to be = 1.6 GeV, = 4.7 GeV and Mf = 
175 GeV. (Note that a complete prediction of the top quark cross section would require 
the incorporation of the axial contribution.) In order to study the sensitivity of the results 
to the renormalization scale, /x^ has been chosen as Mq [Q = c,b,t] dashed curves), 
4Mq (solid curves) and s (dotted curves). The results are relatively stable against this 
variation. This agreement, despite the appearance of large logarithms In s/Mq for /i^ = 
Mq, is a consequence of the NNLO approximation valid at the high energy end. The 
largest sensitivity towards the choice of yU^ is observed in the intermediate and low energy 
region for charm production, where ag is large and the corrections are enhanced by large 
contributions proportional to In f3 and 1//3. For comparison Rc, Rb and Rt are also plotted 
in the Born (wide dots) and leading order approximation (dash-dotted, /i^ = {2MqY). 
The cubic corrections in are rather small so that we do not show the NLO corrections 
separately. For the case of the charm and bottom quark the "remainder" of the Coulomb 
singularity is still visible and both Rc and Rb raise for -^i — > 5 GeV and y/s 11.5 GeV, 
respectively. For the top quark, however, even 10 GeV above the nominal pair production 
threshold a singular behaviour is not visible and Rt decreases as a/s approaches 360 GeV. 

To study the dependence on the input parameters Mq and a^^^(M|), we adopt the 
choice fx^ = ^Mq and vary af\M'^) within the currently quoted range ±0.03 (see 
Fig. ^(a)) and the quark masses in the range as indicated in Fig. ^(b). These changes 
indicate the present uncertainties. For R^ and Rb the higher values of the quark masses 
also lead to larger cross sections — a consequence of the enhanced Coulomb forces for 
fixed v^. For the top quark the phase space effect is still dominant and the cross section 
for the choice Mt = 180 GeV is smaller than for Mt = 170 GeV. 
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1.494 


1.465 


1.447 


1.436 


1.429 


1.425 


0.118 


1.6 


R(s) 


3.696 


3.629 


3.594 


3.572 


3.557 


3.546 


3.542 


0.118 


1.8 


Rc{s) 


1.598 


1.519 


1.480 


1.458 


1.444 


1.434 


1.430 


0.118 


1.8 


R{3) 


3.741 


3.654 


3.609 


3.583 


3.565 


3.552 


3.547 


0.121 


1.4 


Rc(s) 


1.527 


1.480 


1.456 


1.442 


1.433 


1.427 


1.424 


0.121 


1.4 


R{s) 


3.677 


3.621 


3.591 


3.572 


3.559 


3.549 


3.545 


0.121 


1.6 


Rc{s) 


1.569 


1.503 


1.471 


1.452 


1.440 


1.432 


1.429 


0.121 


1.6 


Ris) 


3.719 


3.644 


3.605 


3.582 


3.566 


3.554 


3.549 


0.121 


1.8 


Rc{s) 


1.622 


1.531 


1.488 


1.463 


1.448 


1.438 


1.434 


0.121 


1.8 


R(s) 


3.771 


3.672 


3.623 


3.593 


3.574 


3.561 


3.555 



Table 1: Numerical values for R{s) and Rc{s). The scale fi"^ = s has been 
chosen. For the evaluation of Rc{s) Eq. (|26| ) has been used. R{s) is the proper 
sum of Rudsis) (Eq. (D) and R^s) with = 4.7 GeV. 



It is now possible to give a prediction for R{s) in the energy range above ^/s ^ 2 GeV 
with the exception of small windows of 1 to 2 GeV above the thresholds for open charm 
and bottom production, respectively. Below the charm meson threshold at ^/s = 2M/) 
the charm quark is treated as heavy and Eqs. (0) and are used with rii = 3, and with 
a^'^^ and replaced by a^^^ and respectively. For 5 GeV <, ^/s <, 10.5 GeV Eqs. (||) 
and pip are directly applicable with nj = rii + 1 = 4. For the case with external charm 
and internal bottom quark an expression analogue to Eq. (^) is used. Above the bottom 
threshold the sum over i in Eqs. (|I|) and includes also the charm contribution and 
consequently = 4 has to be chosen, a^^^ and have to be replaced by af''^ and M^. 
The dominant charm mass terms are included through the leading terms in the M^/s 
approximation. The quadratic charm mass corrections of order are given in Eq. (ESf). 



The corresponding terms up to order read as follows [|19 
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TC 
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V* 


11.5 


12 


13 


14 


15 


40 




(GcV) 


(GpV) 














0.115 


AA 


Rb{s) 


0.372 


0.371 


0.368 


0.366 


0.364 


0.349 


0.115 


4.4 


R{s) 


3.899 


3.895 


3.888 


3.881 


3.876 


3.827 


U. i 10 


A 7 




U.O / z 


U.O / ± 


u.ooy 


yj.'ju 1 


U.ODO 




0.115 


4.7 


R(s) 


3.899 


3.895 


3.889 


3.882 


3.877 


3.827 


0.115 


5.0 




0.370 


0.370 


0.369 


0.368 


0.366 


0.350 


0.115 


5.0 


R(s) 


3.897 


3.894 


3.889 


3.883 


3.878 


3.827 


0.118 


4.4 


Rb(s) 


0.376 


0.374 


0.371 


0.368 


0.365 


0.350 


0.118 


4.4 


R{s) 


3.911 


3.906 


3.897 


3.890 


3.884 


3.832 


0.118 


4.7 


Rb(s) 


0.377 


0.375 


0.372 


0.370 


0.367 


0.350 


0.118 


4.7 


Ris) 


3.911 


3.907 


3.899 


3.891 


3.885 


3.832 


0.118 


5.0 


Rb{s) 


0.377 


0.376 


0.373 


0.371 


0.368 


0.350 


0.118 


5.0 


Ris) 


3.911 


3.907 


3.900 


3.893 


3.887 


3.832 


0.121 


4.4 


Rb(s) 


0.381 


0.378 


0.374 


0.370 


0.367 


0.350 


0.121 


4.4 


Ris) 


3.922 


3.917 


3.907 


3.899 


3.892 


3.837 


0.121 


4.7 


Rb{s) 


0.383 


0.380 


0.376 


0.372 


0.369 


0.350 


0.121 


4.7 


Ris) 


3.924 


3.919 


3.909 


3.901 


3.894 


3.837 


0.121 


5.0 


Rbis) 


0.384 


0.382 


0.377 


0.374 


0.371 


0.351 


0.121 


5.0 


Ris) 


3.926 


3.920 


3.911 


3.903 


3.896 


3.837 



Table 2: Numerical values for R{s) and Rb{s). The scale fi"^ = s has been 
chosen. For the evaluation of Rb{s) Eq. (|26|) with obvious modifications has 
been used. R{s) is the proper sum of Rudsc{s) (obtained from Eq. (^)) where 
in addition charm mass effects of order M^/s have been included (Eqs. (pTf ) 
and (0)) and Rbis) (M^ = 1.6 GeV). 





(GeV) 


5 


6 


7 


8 


9 


9.98 


10.52 




Rc{s) 


1.514 


1.437 


1.409 


1.400 


1.397 


1.398 


1.399 


2Mc 


Rc{s) 


1.553 


1.484 


1.450 


1.432 


1.421 


1.414 


1.411 




Rc{s) 


1.553 


1.494 


1.465 


1.447 


1.436 


1.429 


1.425 


2v^ 


Rc{s) 


1.541 


1.493 


1.466 


1.449 


1.438 


1.430 


1.427 



Table 3: Numerical values for Rc{s) for different choices of ^/s and fi. The 
values )(M|) = 0.118, = 1.6 GeV and Mb = 4.7 GeV have been chosen. 



In addition the term 



-7.877 + 0.350n 



f) 



{2i 



which comes from the expansion in M^/ s of diagrams with internal charm loops in hh 
production |l^ has to be added to Rb lis- 
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(GeV) 


11.5 


12 


13 


14 


15 


40 


Aft 


Rbis) 


0.386 


0.381 


0.375 


0.370 


0.366 


0.349 


2Mb 


Rb{s) 


0.379 


0.377 


0.373 


0.370 


0.367 


0.349 


V~s 


Rtis) 


0.377 


0.375 


0.372 


0.370 


0.367 


0.350 


2v^ 


Rb{s) 


0.370 


0.370 


0.369 


0.368 


0.366 


0.350 



Table 4: Numerical values for Rb{s) for different choices of ^/s and fi. The 
values H^l) = 0.118, = 1.6 GeV and Mb = 4.7 GeV have been chosen. 



In Fig. R{s) is plotted versus y/s for the three energy intervals. Solid and dashed 
curves show our prediction from perturbative QCD, where the extreme values from the 
variation of a^^''(M|) and the masses are considered and the scale jj? = s has been 
adopted. For comparison we also plot a recent compilation of the available experimental 
data [^] from the inclusive measurements of R. The central dotted curves correspond 
to the mean values, upper and lower dotted curves indicate the combined statistical and 
systematical errors. 

In Tabs. |1] and |^ numerical values are listed for the energy range between the charm 
and bottom threshold and above the bottom threshold, respectively. For completeness the 
terms of order al to -Rught are also taken into account. They are obtained from Eq. (|25| ) 
neglecting the mass corrections with riy = 4 (Tab. ^ and nj = 5 (Tab. |^), respectively. 
In Tabs. |^ and ^ the dependence of Rc and Rh on the renormalization scale /i is displayed 
with Mq and fixed to their central values. 

It is evident that the uncertainties in the prediction are far below the experimental 
errors in this "low" energy region. These results could therefore be used to fit the cur- 
rently available data and to allow for an improved input into the analysis of the running 
electromagnetic coupling constant a. 



4 Hadronic vacuum polarization and the running 
electroweak coupling 

One of the important ingredients of electroweak precision tests is the effect of the hadronic 
vacuum polarization on the running of the electromagnetic coupling. Using a dispersion 
relation, it is expressed ||21| through R{s) 



A<,(Mi) = -— ^Re / ds- ' ' . (29) 

i-K Jiml S[S — M| — It) 
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0.115 


0.118 


0.121 


Energy r9.iigG 




















2 on — 2 'sn OoV 




7.54 






7.58 






7.62 








6.12 






6.15 






6.18 




3.00 - 3.73 GcV 




"7 OV 
i .li 






■7 Qfl 

( .oU 






/ .OO 






1.4 


1.6 


1.8 


1.4 


1.6 


1.8 


1.4 


1.6 


1.8 


Energy rSyiige 




















3.73 - 5.00 GeV 


9.73 


9.74 


9.74 


9.77 


9.77 


9.77 


9.80 


9.80 


9.80 


(without cc) 




















5.00 — 5.50 GeV 


5.37 


5.41 


5.46 


5.40 


5.44 


5.50 


5.42 


5.47 


5.54 


5.50 — 6.00 GeV 


4.88 


4.91 


4.94 


4.89 


4.93 


4.97 


4.91 


4.95 


4.99 


6.00 — 9.46 GeV 


25.30 


25.37 


25.46 


25.37 


25.45 


25.55 


25.44 


25.53 


25.63 


9.46 — 10.52 GeV 


5.88 


5.89 


5.90 


5.90 


5.90 


5.91 


5.91 


5.92 


5.93 


5.50 — 10.52 GeV 


36.06 


36.17 


36.30 


36.16 


36.28 


36.43 


36.26 


36.39 


36.55 


A/fc (GeV) 


4.4 


4.7 


5.0 


4.4 


4.7 


5.0 


4.4 


4.7 


5.0 


Energy range 




















10.52 - 11.50 GeV 


4.94 


4.94 


4.95 


4.95 


4.95 


4.96 


4.96 


4.96 


4.97 


(without 6b) 




















11.50 - 12.00 GeV 


2.61 


2.61 


2.61 


2.62 


2.62 


2.62 


2.63 


2.63 


2.63 


12.00 - 12.50 GeV 


2.51 


2.51 


2.51 


2.51 


2.51 


2.51 


2.52 


2.52 


2.52 


12.50 - 13.00 GeV 


2.41 


2.41 


2.41 


2.42 


2.42 


2.42 


2.42 


2.42 


2.42 


13.00 - 40.00 GeV 


72.78 


72.80 


72.81 


72.90 


72.92 


72.94 


73.03 


73.04 


73.06 


40.00 - oo GeV 


42.61 


42.61 


42.62 


42.67 


42.67 


42.67 


42.73 


42.73 


42.73 


12.00 - oo GeV 


120.31 


120.33 


120.34 


120.50 


120.52 


120.54 


120.70 


120.72 


120.74 



Table 5: Contributions to I\q\2,^{M'^ x 10^ for different energy regions, 
quark masses and q;^^^(M|). The scale = s has been adopted. Below the 
charm threshold the value = 1.6 GeV is chosen. In the energy ranges 
3.73 - 5.00 GeV and 10.52 - 11.50 GeV for R{s) the formulae vahd below the 
corresponding quark threshold have been used. 



and contributes together with the well known leptonic contributions to the running of the 
electromagnetic coupling 



a(0) 

l-^"had(s) - Aaiep(s) 



«(^) = : . (5)// . 7-, (30) 



where a{0) = 1/137.0359895 is the fine structure constant. Top quark contributions are 
not considered in this context. Also QED corrections are not included as they are of the 
order of a few per mill only. 

The detailed phenomenological analyses in Refs. p2| , p3| have made use of the full 



set of data obtained by many different experiments for energies from just above the 
two pion threshold up to 40 GeV. Although different prescriptions for the interpolation 
have been used in the different papers, the most recent results are in fair agreement. 
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Mc 


2Mc 




2y^ 


Encrg 


y range 










5.00- 


- 5.50 GeV 


O.O 1 


5 44 


5 44 




5.50- 


- 6.00 GeV 




4 92 




4 92 


6.00- 


- 9.46 GeV 


25.10 


25.35 


25.45 


25.46 


9.46 - 


- 10.52 GeV 


5.85 


5.88 


5.90 


5.91 


5.50- 


- 10.52 GeV 


35.80 


36.14 


36.28 


36.28 


M 


Mb 


2Mt 




2v^ 


Energ 


y range 










11.50 


- 12.00 GeV 


2.62 


2.62 


2.62 


2.62 


12.00 


- 12.50 GeV 


2.52 


2.51 


2.51 


2.51 


12.50 


- 13.00 GeV 


2.42 


2.42 


2.42 


2.41 


13.00 


- 40.00 GeV 


72.89 


72.91 


72.92 


72.92 


40.00 


- oo GeV 


42.63 


42.65 


42.67 


42.67 


12.00 


- oo GeV 


120.45 


120.49 


120.52 


120.52 



Table 6: Contribution to Aa^^^^{M^) x 10^ for different energy regions and 
different choices of fi for the most heavy quark contribution (for the remaining 
hght quarks the scale jj? = s has been adopted). The values af^M"^) = 0.118, 
Mc = 1.6 GeV and Mf, = 4.7 GeV have been chosen. 



Energy range 


this work 


Ref. y 


5.00 - 9.46 GeV 
12.00 - 40.00 GeV 


(35.32 - 35.81) ± 0.4 
(77.82 - 77.84) ± 0.2 


32.63 
79.22 



Table 7: Comparison of Aa{^3^(M|) x 10^ evaluated in this work with the 
data analyzed in |^. The scale fi^ = s has been chosen for the respective light 
quark contributions and varied in the range between Mq and s for the most 
heavy quark. The errors arise from the variation of M^, and a^^-*(M|) (see 
Tab. I). 



In the high energy region (typically above 40 GeV p2|) the prediction for R{s) based 
on perturbative QCD with massless quarks has been employed. A significant part of 
the final error originates from the region where perturbation theory should be reasonably 
valid: the light quark continuum, say above 2 GeV, the continuum above the charmonium 
resonances and below the bottom threshold and the region between 12 and 40 GeV, i.e. 
above the Upsilon resonances. In view of the results presented in the previous chapters 
one may employ perturbative QCD also in these regions. This might lead to a reduction 
of the error, albeit at the price of a more pronounced dependence on perturbative QCD. 
(For early studies along this line see also ||2^, |23[|.) 

In Tab. ^ and ^ the contributions to Aa[j^^^ are displayed for a variety of input pa- 
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rameters and renormalization scales. Since the validity of our perturbative treatment is 
more doubtful just above the respective charm and bottom thresholds, the contributions 
from a variety of intervals are displayed separately. Once improved data are available, 
these theoretical numbers could be replaced by more precise experimental ones. In Tab. ^ 
our results are compared to the analysis of [E^ for two energy ranges. For the lower 



energy interval from 5.00 — 9.46 GeV our QCD prediction for Aaj^^^j is bigger than the 
value obtained by including experimental data by about 9%. In contrast, for the energy 
interval from 12 — 40 GeV, perturbative QCD gives a value slightly smaller than the one 
obtained from integrating the experimental R{s) values. This behaviour is already clear 
from Fig. |^, where we indicate the range of experimental data in comparison to R{s) from 
perturbative QCD. 

If one subtracts the narrow Upsilon resonances, QCD can be applied even up to 
the threshold for B meson production. In fact, recent experimental investigations at 



10.52 GeV |2J] are nicely consistent with perturbative QCD. The additional contributions 
from the continuum cross section in the intervals between 3.73 GeV and 5.00 GeV (without 
cc) and between 9.46 GeV and 10.52 GeV (without bb) are also listed in Tab. |[ For the 
perturbative contributions through the range from ^/s = 3 GeV to oo without charmonium 
and open charm contributions below 5 GeV and T and open bottom contributions below 
11.5 GeV one finds 

Aafadlpert = [(186.27 - 186.87) ± 0.70] x 10"^ (31) 

where the range is due to the variation of fi for the heavy quark contribution and the 
error due to uncertainties in the parameters. A more detailed discussion of the impact of 
these calculations will be given elsewhere. 



5 The region for small (3 — closer to the threshold 

In the high energy limit, say down to s ~ 8Mq, the cross section in 0{as) and 0{al) 
is well described by the massless approximation plus the leading terms of the expansion 
in Mq/s up to {Mq/s)^ [^, ^. The bulk of the large logarithms is resummed by taking 
n"^ = s for the renormalization scale and adopting the MS definition of the running 
mass. The fixed order result as given above is evidently adequate in the intermediate 
energy region, with the requirement that CpiTas/P is not yet too large, i.e. safely away 
from the threshold regime, where the conventional multi-loop expansion breaks down. 
Otherwise some care has to be taken to control the higher order terms proportional to 
[CpT^a/ (3)"' with n > 0. As far as dominant and subdominant contributions of this sort 
are concerned, their structure is understood from nonrelativistic considerations and will 
be briefly outlined in the following for the case Cfttc^s ^ (3 1. 

Let us in a first step discuss those terms which are multiplied by the colour factor Cp 
only and which are relevant for QCD and QED. For the dominant contributions to the 
cross section in the nonrelativistic limit, often called "Sommerfeld factor" in literature, 
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the leading terms in an expansion in xs = CpT^ots/ (3 are given by 



R 



thr 



3 pxs 



2 1 - e-^s 



Xs Bix 



B2X% ^ Bsx% 



TT 2/5 



Cp TTo + C'l' 



4! 

TT 



6! 

4 



12/5^ 



+ . . 



(32) 



where i?„ are the Bernoulh numbers: Bi = 1/6, B2 = 1/30, B^ = 1/42, ... .It should be 
noted that the Sommerfeld factor is entirely of long-distance origin and proportional to 
the imaginary part of the nonrelativistic Green function for the Coulomb potential p6| , 
i.e. governed by the continuum Coulomb wave function. These terms are predicted from 
a consideration, where the QQ production process is decomposed into a short distance 
part (to be eventually corrected by short distance QCD corrections) and a long distance 
part, which is governed by the Coulomb wave function, in other words, by the imaginary 
part of the nonrelativistic Green function for the Coulomb potential |]27| , ^ . 

Resummed and fixed order results have to coincide in the region of small xs- Thus it 
is instructive to compare the Sommerfeld factor and the sum of the first three terms in 
the expansion (^) (corresponding to Born, one- and two-loop contributions, respectively) 
for various values of X5 which are not too much larger than one. It is remarkable that the 
sum of the first three terms in Eq. (|32|) provides an excellent approximation not only for 
small values of xs, but even up to X5 = 2 with a relative deviation of less than 1%. Even 
for Xs = vr, corresponding to Cpas/P = 1, the deviation amounts to 3.3% only. 

In addition, based on the consideration that the cross section for nonrelativistic ener- 
gies can be decomposed into long- and short- distance contributions one obtains in QED 
an additional correction factor which comes from relativistic momenta involved in the 
transverse photon exchange. This factor, which is quite familiar from the single-photon 
annihilation contributions to the positronium hyperfine splitting |2^ and from the correc- 
tions to quarkonium annihilation through a virtual photon can be derived from the 
threshold behaviour of the one- loop corrections [|l^ . In combination one thus anticipates 
the following behaviour 



R 



thr 



3QlP 



3-/32 




1 -4C, 



TT 



(33) 

The inclusion of the subleading (3'^ terms in the combinations (3 — (3^) and (1 + (3"^)/ (3 is 
suggested from the structure of Eq. (P) for small (3. An explicit proof for this structure at 
the NLO level can be found in |^T|. Their interpretation as long distance contributions is 
furthermore motivated by the appearance of the logarithmic terms with the same structure 
in the CpTrii contributions listed below (see |]15[ and Eq. (40) of [§]). Even the impact 
of the running of the coupling constant can be included in this nonrelativistic line of 
reasoning. In QCD the running of ag is induced by terms proportional to Tn^ and Ca- 
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For a prediction of the hadronic cross section very close to threshold, i.e. in the regime 
/3 < CpTTasi the factorization into short- and long-distance contributions analogous to 
Eq. (|33|) is desirable, incorporating also the new information from the NLO calculation. 
Such an analysis requires the resummation of long-distance contributions to all orders and 
shall not be carried out here^. However, the leading terms of the perturbative series for 
small (3 provide the basis for such a resummation and will be collected in the following. 

Let us now recapitulate the threshold behaviour of the various ingredients for Rf.: 



4°) 
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2 ' 
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12 

TT 
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+ 



(34) 
(35) 



,7r2(l + /32) 
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/5 
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TT 



M, 



Ml 
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+ 



4, 

+ 3^^m! + 



11 



+ . . . 



44 8 



^ln-^-,ln^ + ^-,C(2) 



9 



+ 



(36) 
(37) 
(38) 
(39) 



Terms of order (3"^ (modulo In /?) are neglected. The order /3 terms in Eq. (^) is taken 



from |3^. Subleading terms in Eq. (0), symbohzed by c, have not been calculated yet 
analytically. Assuming a linear dependence on /?, an estimate for the constant c can be 
extracted from the numerical analysis in |p. Considering the expansion of 11 different 
Fade approximations which show a quite stable behaviour near threshold one obtains 



24 ±5, 



(40) 



where the error is estimated by taking into account the variation of the predictions for the 
constant c from the different Fade approximations. This numerical result for c is fairly 
large, in particular when compared to the corresponding constant 11/9 in the CpTrii 
term, Eq. (H). 

In total one thus gets|: 



Rr 



Ql 



TT 



^ A resummation of this sort for the QED contributions, i.e. inchiding also the order f3 terms in 
Eq. (B3), can be found in p3]. 

^ The contributions from SR^j^ are suppressed by the power /M^, not Coulomb enhanced and thus 
ignored in the following. 
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(41) 



The renormalization group invariance of this result is apparent: the dependence of i?7VA, 
i?cg, and Rcc properly compensates the n dependence of the 0{as) term. As stated above, 
the energies considered in this section will be of order Mc, implying that ln(s/M^) is not 

a large quantity. The transition from aP^^ to af'^ is thus legitimate and easily achieved 
by absorbing the last term of Eq. (|11|) in the order as expression. 

The logarithmic singularities and the constants in Rna and Rcq which are leading in 
(3 can be absorbed in the terms of order if the MS coupling constant is replaced by the 
coupling governing the potential 



VqcD^r) 
ay(g^) 



-AtiCf 



a 




(42) 



with = (3'^s. This is apparent once the sum of the Born cross section plus higher order 
corrections is rewritten as follows (with xy = CpT^Oivi^P'^s) (1 + /J^)//?): 



(3) 



Rc = i + ^-4C^" 

Z n 



12 ^2 77 




2 r 



3 I 4 12/ 4 ' 



11 



+ . . . 



(43) 



In the transition from Eq. (^) to (|43| ) we have freely dropped terms of order a^. It 
is evident that the scale in the correction term from hard transverse gluon exchange 
proportional to AC pOi^^^fi^) / is of order Mc, with. fj.BLM = e" 

-11/24 j^^^ ^ 0.63Mc suggested 



by the BLM prescription However, as noted before, the corresponding constant 

c ~ 24 in the non-abelian term is markedly different from what would be expected from 
the BLM procedure. 

As stated above these results are strictly applicable in the limit nCpas ^ /3 <C 1 only. 
Nevertheless, Eq. ( |^ ) provides an important input for the determination of the cross 
section very close to threshold and even for bound state energies, i.e. for \P\ < CpTcas, 
because it contains all 0{a^) short- distance effects relevant for the nonrelativistic regime. 
These short-distance effects, which are specific for the single-photon annihilation process 
involving massive quark- ant iquark pairs, are universal for |/?| <^ 1 regardless whether \(3\ 
is smaller or larger than Cpnas- In this respect Eq. ( ^3]) even provides an important 
result for the investigation of leptonic decay widths of the and the T families, and for 
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QCD sum rules for the bb system. A discussion of the latter subjects, however, is beyond 
the scope of this paper. 

6 Summary and Conclusions 

Predictions for the cross section of massive quark production in e^e~ annihilation are 
presented which are accurate to order a^. Their range of validity extends from high 
energies down close to threshold, i.e. to center of mass energies of about 5 GeV, 11.5 GeV 
and 2Mt + 12 GeV for charm, bottom and top respectively. Inclusion of the leading 
and subleading terms of order proportional to Mq/s allows to connect smoothly the 
NNLO prediction at high energies with the NLO prediction in the intermediate and "low" 
energy range. The NLO corrections are sizable and must be taken into account to achieve 
a prediction with an accuracy better than 10%. The stability of the prediction against 
variations of the renormalization scale and of the input values for the quark masses and 
as has been tested. Even fairly extreme assumptions about the renormalization scale lead 
to moderate variations in the case of charm and to negligible variations for the heavier 
quarks. The same holds true for the dependence on the input parameters, the quark 
masses and the strong coupling constant. Only a change in the charm quark mass has a 
clearly visible effect on the cross section. Within the large experimental errors, theoretical 
and experimental results are well consistent. The theoretical results can now be used to 
perform an universal fit to the data, including the lower energy regions. In view of the 
sparse data with large errors in the region from 2 — 3.73 GeV, from 5 — 10.52 GeV and 
from 11.5 — 40 GeV one could also consider to use the predictions based on perturbative 
QCD to arrive at a more precise value for the running QED coupling at the Z boson 
mass. A significant reduction of the uncertainty could be obtained. 
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